Mott transition between a spin-liquid insulator and a metal in three dimensions 
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We study a bandwidth controlled Mott metal-insulator transition (MIT) between a Fermi liquid 
metal and a quantum spin-liquid insulator at half-filling in three dimensions (3D). Using a slave rotor 
approach, and incorporating gauge field fluctuations, we find a continuous MIT and discuss the finite 
temperature crossovers around this critical point. We show that the specific heat C ~ Tlnln(l/T) 
at the MIT and argue that the electrical transport on the metallic side near the transition should 
exhibit a 'conductivity minimum' as a function of temperature. A possible candidate to test these 
predictions is the 3D spin liquid insulator Na4lr308 which exhibits a pressure-tuned transition into 
a metallic phase. We also present the electron spectral function of Na^^Os at the transition. 



Introduction. — Recent experiments suggest that 
Na^^Os at ambient pressure may be the first exam- 
ple of a three-dimensional quantum spin-liquid Mott- 
insulator[l]. This material has Ir local moments residing 
on the hyper-kagome lattice which is a three-dimensional 
(3D) network formed by corner sharing triangles. Of the 
various theoretical proposals for the spin liquid phase ob- 
served in this material [2-6], the data may be most con- 
sistent with a particular gapless spin liquid whose low 
energy excitations are charge-neutral spin-1/2 fcrmions, 
called spinons, which live on 'spinon Fermi surfaces' in 
momentum space [4, 5] . Remarkably, very recent exper- 
iments show that this material undergoes a transition to 
a metallic state under pressure [7]. This indicates that 
Na4±r308 may be a 3D analog of the well-studied layered 
triangular lattice organic material k— (ET) 2 Cu 2 (CN) 3 [8] 
which is insulating at ambient pressure but undergoes 
a Mott transition to a metal under moderate pressure. 
Theoretically, the proximity to the Mott transition has 
been argued to lead to a 2D spin liquid with a 'spinon 
Fermi surface' in k-(ET) 2 Cu 2 (CN) 3 [9, 10]. 

In this Letter, motivated by the recent experiments 
on Na4iX308, we present a theory of a continuous Mott 
transition between a spin-liquid with a spinon Fermi sur- 
face and metal in 3D. We use a similar theoretical frame- 
work as developed earlier by one of us[ll] for a 2D Mott 
transition between a Fermi liquid and a spin liquid Mott 
insulator in relation to k— (ET) 2 Cu 2 (CN) 3 . In 2D, sev- 
eral interesting results such as a universal jump in the 
residual resistivity and the emergence of marginal Fermi 
liquids in an intermediate crossover scale were found. In 
this work, we show that a continuous Mott transition is 
indeed possible in 3D and moreover it is characterized by 
weaker singularities in various physical quantities when 
compared to the 2D case. In contrast to 2D, right at the 
Mott critical point in 3D, the system is insulating. In- 
deed we argue that the electrical conductivity a(T) in the 
vicinity of the Mott quantum critical point (QCP) should 
exhibit the following behavior: (i) In the 'quantum crit- 
ical' (QC) regime, T > T*, we expect a(T) ~ T (with 
log corrections). Here the crossover scale T* ~ \g — g c \i 
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FIG. 1: (a) Schematic phase diagram of the 3D metal insula- 
tor transition. Finite temperature crossovers at the temper- 
atures T* and T** , shown by dashed lines, separate regions 
distinguished by the specific heat and conductivity behaviors, 
as discussed in the text ( "C-SL" labels the Charge fluctuation 
renormalized Spin Liquid), (b) Resistivity and (c) conduc- 
tivity along the lines - A, B, and C - in panel (a). 

where g is the parameter used to tune the transition (for 
instance the pressure, or the strength of the Coulomb 
repulsion U/t). g c is the value of the parameter at the 
zero temperature Mott QCP. (ii) On the insulating side, 
for T < T* , the conductivity is thermally activated, (hi) 
Finally, on the metallic side, for T ~ T* , we argue that 
cr(T) will have a 'conductivity minimum' before it grows 
and saturates at low temperature. Equivalently there will 
be a resistivity peak at a non-zero temperature. This is 
sketched in Fig. 1(b) and (c). This resistivity peak is the 
most striking result of this paper. 

As in the 2D case the finite temperature crossovers 
near the Mott transition are characterized by multi- 
ple energy scales T* and T** as depicted schematically 
in Fig. 1(a). As noted above the scale T* determines 
the crossover in transport properties associated with the 
Mott transition. However there are only weak changes in 
physical properties across the scale T** ~ \g — g c \* . For 
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instance for T 3> T**, in the regimes labelled quantum 
critical (QC), 'anomalous metal' and 'charge fluctuation 
renormalized spin liquid' (C-SL), the specific heat be- 
haves as C ~ Tlnln(f/T). For T <C T** in the Mott 
insulator C ~ Tln(l/T) while in the same regime on the 
metallic side we recover the Fermi liquid result C ~ T. 
With an eye toward application to Na^^Os, we also 
present predictions for the T = electron spectral func- 
tion and tunneling density of states at the Mott transition 
for a tight binding model appropriate to this material. 

Microscopic Model and Field Theory. — To understand 
the universal features of this kind of 3D Mott transi- 
tion we start with a single band Hubbard model at half- 
filling on a non-bipartite lattice. H = —tJ2(ij) Aa c j° + 
T Si n i~l J Si n «' where c\ a is the creation operator for 
the electron at site % with spin cr, and rij = c\ a Ci„ is 
the electron number operator (repeated spin indices are 
summed over). This model Hamiltonian is expected to 
have a metallic Fermi-liquid ground state for small on- 
site repulsion U/t <C 1 , and to have a Mott insulating 
ground state for U/t 3> 1. (Henceforth, we set t = 1.) 
The arguments of Ref. 9 suggest that a spin liquid state 
with a spinon Fermi surface is a likely ground state just 
on the insulating side of Mott transition. 

The Mott transition between a Fermi-liquid metal and 
such a spin-liquid insulator is most conveniently accessed 
in a slave rotor formalism [12], in which the electron op- 
erator is written as a product, c\ a = e Wi f} a , of a charged 
spinlcss 'rotor' field <fii = and a spin-1/2 charge- 
neutral 'spinon' fi a . In order to eliminate unphysical 
states from the enlarged spinon-rotor Hilbert space, we 
must impose the constraint n\ + n{ = 1 at each site. 
Upon coarse-graining, this theory can be recast in the 
form of a 'gauge theory' which consists of a dynamical 
U(l) gauge field et M = (ao, a) coupled to, both, the spinon 
and rotor fields [10]. We begin by writing down a contin- 
uum theory to describe this physics in terms of a complex 
bosonic field <fi ~ e lS , a fermionic field tjja, and a gauge 
field (fy, with action S — dr J d 3 r C, where 

C = Cb + Cf + Cg + Cbf 

C b = |(a /1 + ia„M 2 + m 2 M 2 + uM 4 

t 1 2 

Cf = ipl(d T - ia - Hf)tp<T + |(V-ia)?/v 

2m f 

Cbf = A|^| 2 |0| 2 (1) 

These terms are all allowed by the symmetries of the sys- 
tem, and thus will generally arise in the coarse-graining 
process. Note that the boson action is particle-hole sym- 
metric, as appropriate for a coarse-grained version of the 
slave rotor theory at half-filling. 

At mean field level, we can ignore fluctuations of the 
gauge field a^. In the regime m 2 > 0, the charged rotor 
field is gapped, so that {4>i) = 0, while the spinons form 



a gapless Fermi surface. This is the spin liquid insula- 
tor. On the other hand, the metallic phase at to 2 < 
corresponds to a Bose-Einstein condensate of rotors with 
(<f>i) ^ 0. The spinon Fermi surface now becomes an 
electron Fermi surface with quasi-particle weight propor- 
tional to \(4>i)\ 2 ■ Therefore, the transition to the metal is 
tuned by the boson gap parameter m 2 , with m 2 oc U—U c . 
The MIT, in this mean field theory, is thus simply a Bose 
condensation of rotors with the spinon Fermi surface be- 
ing continuously connected to the electron Fermi surface, 
and hence it is a continuous transition. 

We now consider a full theory of this Mott transition 
obtained by analysing fluctuation fluctuation corrections 
to all low-energy degrees of freedom (spinons, rotors and 
gauge fields) in a self-consistent fashion. As the methods 
have already been elaborated in the 2D case in Ref. 11 
we will be brief and mainly state the results. 

We begin by considering the boson interaction term 
u\4>\ together with the fermion-boson density interaction 
Cbf, in the absence of the gauge field a M . (The gauge field 
is reinstated below). Upon integrating out fermions (see 
Fig. 2(e)), this yields an effective boson-boson interaction 

k, k' , q 

where u = u + Nq\ 2 it f(q,iuj n ), No is the fermionic den- 
sity of states, and 7r/(g, iui n ) is the density-density polar- 
ization function of the fermions. For \q\ <C 2/c/, and 
in the limit |w n | <C Vjq, 717 — > 2 — ^j^- Under a 
naive rescaling of space and time with dynamical crit- 
ical exponent z = 1, £ is marginal by power count- 
ing at the critical fixed point of the bosonic theory, just 
like the conventional u|0| 4 term. The fate of this term 
can be determined by a standard perturbative renormal- 
ization group analysis. Equivalently we generalize the 
theory to Nj, boson flavors, and consider the Nj, = 00 
limit of the theory. In this limit, all diagrams can be 
summed exactly to yield a renormalized bosonic interac- 
tion Uren(<7> w ) = u^ 1 (q 7 uj) +ln(A / ^uj 2 + q 2 + m 2 ). Near 
the QCP, the second term dominates, and u Ton — ► 0. 
Hence, C is (dangerously) marginally irrelevant. 

Now we reinstate the gauge field and first consider its 
renormalization from integrating out the matter fields. 
We work in the Coulomb gauge V • a — 0, in which ao 
and a decouple in C g . Then, the longitudinal field ao 
is screened by the gapless fermions, and therefore can 
be ignored. On the other hand, the transverse part 
of a remains unscreened, but it is strongly renormal- 
ized by the matter fields, whose contribution dominates 
over the bare gauge field action. In the Random Phase 
Approximation (RPA), the renormalized inverse trans- 
verse gauge field propagator is given by D~ (q, iv n ) = 
nj(q, iu n ) +IIh(q, iv n ) where and lib are fermion and 
boson polarization functions. 
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FIG. 2: Dominant fluctuations. The solid and dashed ar- 
rows represent fermions and bosons, respectively, (a) and (b) 
are the gauge polarization functions arising from integrating 
out fermions and bosons, respectively, (c) Spinon self-energy, 
where the double wavy line is the dressed gauge propagator, 
(d) Boson self-energy, (e) Residual bosonic interactions. 

The fermionic contribution, shown in Fig. 2(a), is ILf = 
7o + Xoq 2 '■ The term proportional to 70 describes the 
Landau damping due to the gapless spinons and xo is 
the diamagnctic susceptibility of the spinons. This form 
can be further justified in the large Nf limit, where Nf 
is the number of flavors of the fermions. The correc- 
tions appearing in a 1/Nf expansion do not change this 
form[13]. The bosonic contribution shown in Fig. 2(b) for 

m 2 > (or U > U c ) is II 6 = ^ In (g2+ "° 2)1/2 , where m 

2 

is a non- universal constant, leading to lib = Ss hi — 

2 

at the critical point U = U c and Ilf, « ^ n 7? f° r 
U > U c . On the other hand, when m 2 < or U < U c , 
II5 w p s where p s oc |{</>)| 2 is the 'superfluid stiffness' 
of the bosons. Thus the bosons determine the different 
forms of the gauge field propagator in various regimes. 

It is readily seen that the boson self-energy £&(q, u>), 
given by Fig. 2(d), acquires only analytic corrections in 
q and u and thus the bosons are not renormalized in an 
essential way by the gauge fluctuations. The arguments 
in Ref. 11 now imply that the boson sector is basically 
not affected by the gauge fluctuations in the scaling limit 
in the vicinity of the critical point. 

The fermion self-energy shown in Fig. 2(c) can be 
computed by using the gauge field propagators obtained 
above. Near the spinon Fermi surface, |k| s=s fc/, this 
leads to £j ~ u In In l/\u>\ ln ^jnj\ at the critical point 
U = U c and £y ~ u>lnl/|w| + in the spin liquid 

phase for U > U c . 

Finally direct calculation shows that the c a — > 4> + f a 
vertex is not singularly enhanced at low frequency. Thus 
we may calculate the singular part of the physical elec- 
tron Greens function by simply convolving the <ft and 
f a Greens functions. In the metallic phase, we get 
the electron Green's function Q e = \{(j))\ 2 Qf = u _f" +iS 



With 6 *P ^ MR and Z e - lnln'l/p. - Usin S P* K 

(U c — U)ln u _ v , the effective mass of the electrons 
m*/m e ~ lnln 1/(U C — U) diverges very weakly. The 

• 1 rr (U c -U) ha. tj-^jt ,. ... 

quasi-particle weight Z e ~ lnln diminishes as 

one approaches the QCP. 

Thus we see that there are only weak calculable loga- 
rithmic corrections due to fluctuations on top of the slave 
rotor mean field theory of the Mott transition in 3D. This 
enables us to make a number of firm predictions about 
the universal properties of the transition. 

Scaling, Finite Temperature Phase Diagram, and Spe- 
cific Heat. — As in 2D the finite T crossovers will involve 
multiple energy scales originating in the different space- 
time scaling of the gauge field and the bosons. Note that 
the mean-field transition of bosons at finite temperature, 
below which p s ^ 0, becomes a crossover in the full the- 
ory and represents the onset of the charge coherence at 
this temperature scale. 

The behavior of the charge excitations is governed by 
T ~ io q scaling [z = 1 theory). In either phase for 
T above a crossover scale T* ~ \U — U c \ 1,/2 , the charge 
bosons arc in their quantum critical regime. In the in- 
sulator the scale T* may also be identified with the zero 
temperature charge gap. 

The scaling behavior of the coupled spinon-gauge sys- 
tem is, however, determined by T ~ u> ~ q 3 scaling 
(2 = 3 theory) up to logarithmic corrections (in principle, 
one should use u> ~ q 3 ln 1/q scaling) . The spinon-gauge 
system thus emerges out of the quantum critical regime 
only at a scale T** ~ \U — U c \^. This crossover scale is 
in principle visible in the specific heat which behaves as 

(Tlnlnl/T T >\U — U c \ 3 ^ 2 
C^ljiTlnl/T T<(U-U C ) 3 / 2 (2) 

[ 72 T T < {U c - U) 3 ' 2 

where 71 oc 1/ ln ryz^j ) and 72 ~ ln ln ^ 3^ . Note 
however that there only are rather weak changes in prop- 
erties across this second crossover temperature T** . 

The crossover scale T* is, on the other hand, visible 
in the electrical conductivity which is readily measured 
in experiments. It can be obtained from the Ioffe-Larkin 
rule, which states that in a slave-particle formalism, the 
electrical resistivity of the system is the sum of the resis- 
tivities of the fermions and the bosons, i.e. p = pf + pi,. 
For the fermions, wc assume that the conductivity is 
dominated by elastic scattering due to impurities in the 
system, a f = <Jffi. The bosonic conductivity, on the 
other hand, is strongly temperature and pressure depen- 
dent. Inside the Mott insulator, when the bosons are 
gapped (T < T*), the bosonic conductivity is activated 
CTb(T) ~ e~ T l T . On the other hand, inside the metal, 
the bosons are condensed, and ct&(T — > 0) — » 00. In 
the quantum critical regime, standard scaling arguments 
give (Tb(T) ~ Tin 2 y (the log corrections are due to the 
marginally irrelevant boson-boson interactions). As this 
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goes to zero for small T, the boson contribution domi- 
nates over that of the fermions due to the Ioffe-Larkin 
rule. Hence, the critical point itself is insulating. 

These results for the conductivity are summarized in 
Fig. 1(b) and (c). Note that in the metallic phase near 
the critical point, a is not monotonic with temperature. 
At high T, the system is in the (insulating) quantum 
critical fan, and conductivity decreases with decreasing 
T. It is only at low temperatures that the system real- 
izes it is inside the metal, and that conductivity begins to 
rise. Hence, it is possible for a system to be metallic, and 
yet display "insulating" behavior at high temperatures. 
The non-monotonicity becomes more marked for cleaner 
samples, i.e. when o~f_o is large. The experimental ob- 
servation of such a conductivity would provide dramatic 
evidence for the 3D Mott-metal transition discussed here. 

Mean-Field Theory for Hyper-Kagome Lattice. — We 
have shown that the bosons and spinons are only weakly 
affected by the gauge fluctuations near the quantum crit- 
ical point. This suggests that the mean-field theory can 
capture some essential features of the continuous Mott 
transition in 3D. Encouraged by this observation, we 
specialize to the hypcrkagome lattice (of Na^^Os) and 
closely investigate the mean-field theory by approximat- 
ing the ground and excited states via direct products of 
spinon and rotor wave functions |^) = [14]. In 

the mean field analysis [15], we indeed find a continuous 
quantum phase transition at U c ss 6 from a spin-liquid 
insulator to a metallic phase. The electron spectral func- 
tion is given by, 

AeM (3) 
q.a *^ 

x [(i + n(ao -/(£_,))% -£_ q -n q ) 

+ (n(fi q ) + - £k- q + fi q )] 

where is the dispersion of the spinon band a = 1 ... 12, 
fiq is the boson dispersion, and n and / are the Bose 
and Fermi distribution functions, respectively. Here, the 
'form factor' M a (k — q) arises from the spinon wave func- 
tion in band a. At the QCP, we assume the boson dis- 
persion to be given by fi q = 4i^(sin 2 ((fc/2)-|-sin 2 ((7j / /2)-|- 
sm 2 (q z /2)), so that f2 q w «{,|q| for small momenta. For 
illustrative purposes, we choose = 0.5. Deep in the 
metallic phase at T = 0, where the bosons are fully con- 
densed, Ai mot) (k, lo) = J2 a M a (k)S(cu - 

Figure 3 shows the T = spectral function along some 
high symmetry directions of the Brillouin zone (BZ). 
Deep in the metallic phase, the spectral function clearly 
shows twelve bands (after accounting for the four-fold 
degeneracy of the flat band at high energy). Momentum 
variations of the 'form factor' M Q (k) lead to changes in 
the intensity for the various bands as one traverses the 
BZ. In particular, only the lowest band has nonzero in- 
tensity at the T-point. As one approaches the QCP, the 
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FIG. 3: Intensity plot of the electron spectral function along 
high symmetry directions in the metallic phase (left) and at 
the critical point (right). Here F — (0,0,0), R = (n,n,n), 
and M = (ir, n, 0). uj — marks the chemical potential. 

overall bandwidth gets suppressed due to the Hubbard 
U. In contrast to the metallic phase, the spectra at the 
QCP are considerably broadened due to rotor fluctua- 
tions. Finally, as evident from Fig. 3, there is significant 
suppression of spectral weight around the chemical po- 
tential. This 'pscudogap' is also reflected in the tunneling 
density of states obtained by integrating ^4.( cnt ) (k, oj) over 
all k. We find iV t ^"^(w) ~ lo 2 at low frequencies. 

Discussion. — In Ref. 4, it was shown that the specific 
heat of the insulating paramagnetic phase of Na^^Os 
at not-too-low temperatures is well described by a renor- 
malized mean-field theory of a gapless spin liquid with a 
spinon Fermi surface. Such success, however, may seem 
at odds with the fact that the gauge field fluctuations 
would give rise to C /T ~ In 1 /T in the spin liquid regime 
while the renormalized mean-field theory does not neces- 
sarily capture the gauge fluctuation effects at finite tem- 
peratures. In light of the present work, it is possible that 
the spin liquid phase of Na4lr 3 0s may already be near 
the QCP so that the temperature range where the bona 
fide spin liquid phase can be seen is quite narrow. In 
this case, the finite temperature behavior of the specific 
heat may be dominated by the C-SL regime, leading to 
C/T ~ lnlnl/T. This weak singularity would indeed be 
consistent with the apparent success of the renormalized 
mean-field theory. The recent discovery of the insulator- 
metal transition in Na^^Os upon increasing pressure 
[7] may also be consistent with the picture that the spin 
liquid phase discovered earlier [1] is already close to the 
QCP. Furthermore, our analysis suggests that the gap- 
less spin liquid ground state with a spinon Fermi surface 
itself is more likely to occur due to the proximity to the 
Mott transition and the associated charge fluctuations 
while it may face stiff competition from other possible 
phases deep inside the Mott insulator. The confirmation 
of this picture would come from careful comparisons be- 
tween future experiments and the detailed predictions on 
the finite temperature crossover phase diagram, specific 
heat, and conductivity presented in this work. 
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